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We measure the time dependence of the ratio of decay rates for the rare decay D° — > K~^i:~ to 
the Cabibbo- favored decay D^^ K^tv^. Wo use a sig rial of 12.7 x 10''* D° K+tt' decays with 
proper decay times between 0.75 and 10 mean lifetimes. The data sample was recorded with 
the CDF II detector at the Fcrmilab Tcvatron and corresponds to an integrated luminosity of 1.5 
fb~^ for pp collisions at ^/s = 1.96 TeV. We search for _D° — ff' mixing and measure the mixing 
parameters to be Rd = (3.04±0.55) x 10"^, y' = (8.5±7.6) x 10"^, and x'^ = (-0.12±0.35) x 10"^. 
We report Bayesian probability contours in the x'"^ — y plane and find that the data arc inconsistent 
with the no-mixing hypothesis with a probability equivalent to 3.8 Gaussian standard deviations. 
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Since the discovery of the charm quark in 1974 P, [3], 
physicists have been searching for the oscillation of neu- 
tral charm mesons between particle and anti-particle 
states. Such behavior is referred to as "mixing" , as 
first explained in 1955 3] for the meson in terms 
of quantum- mechanical mixed states. Mixing was next 
observed for mesons in 1987 0, i]. The years 2006 
and 2007 have seen landmark new results on mixing: ob- 
servation of Bs mixing from the CDF experiment [6[ and 
evidence for mixing from the Belle [TJ and BABAR Q 
experiments. 

The recent evidence for mixing comes from two dif- 
ferent types of measurements. The Belle Collaboration 
found direct evidence for a longer and shorter lived 
meson, in analogy to the well-known case for K'^ mesons. 
They found significantly different decay time distribu- 
tions for _D" decays to the CP-eigenstates K^K~ and 
7r+7r~ compared to that for the CP-mixed state K^tt^ . 
(In this Letter, reference to a specific decay chain im- 
plicitly includes the charge- conjugate decay.) No other 
experiment has confirmed the evidence for lifetime differ- 
ences among these decays. The evidence for mixing 
found in the BABAR experiment is a difference in decay 
time distribution for K^Tr~ compared to that for 

the Cabibbo-favored (CP) decay D° ^ K-tt+. Such a 
difference depends on the combined effects of differences 
in the masses and lifetimes of the meson weak eigen- 
states. This same measurement was made in the Belle 
experiment Q , but evidence for mixing was not seen. In 
this Letter, we present a new measurement of the same 
D'^ mixing process as used by BABAR for their evidence. 

In the standard model, the decay D'-* A'+tt^ 
proceeds through a doubly Cabibbo-suppressed (DCS) 
"tree" diagram, and may also result from a mixing pro- 
cess (Z?° ^ D'^), if it exists, followed by a CP de- 
cay {D^ K'^TT"). The DCS decay rate depends on 
Cabibbo-Kobayashi-Maskawa quark-mixing matrix ele- 
ments and on the magnitude of SU(3) flavor symmetry 
violation fl3|. Mixing may occur through two distinct 
types of second-order weak processes. In the first, the D'^ 
evolves into a virtual ("long-range") intermediate state 
such as TT+TT^, which subsequently evolves to a D^. The 
magnitude of the amplitude for long-range mi xing has 
been estimated using strong interaction models |11| , but 
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has not been determined using a QCD calculation from 
first principles. The second type of second-order weak 
process is short-range [l^l, with either a "box" or "pen- 
guin" topology. Short-range mixing is negligible in the 
standard model. However, exotic weakly interacting par- 
ticles could enhance the short-ran ge m ixing and provide 
a signature of new physics [l^, H Il5| . 

The ratio R of D° ^ K+tt- to D° K-tt+ decay 
rates can be approximated 1^, 17 1 as a simple quadratic 
function of i/r, where t is the proper decay time and 
T is the mean D'^ lifetime. This form is valid assuming 
CP conservation and small values for the parameters x = 
AM/T and y = AP/2P, where AM is the mass difference 
between the meson weak eigenstates, AP is the decay 
width difference, and P is the average decay width of the 
eigenstates. Under the assumptions stated above, 



Rit/r) 



^y' {t/r) 



y 



{t/r)\ (1) 



where Ru is the squared modulus of the ratio of DCS 
to CP amplitudes. The parameters x' and y' are linear 
combinations of x and y according to the relations 

x' = X cos (5 -f y sin 8 and y' — ~x sin b ^ y cos (5, 

where 8 is the strong interaction phase difference between 
the DCS and CP amplitudes. In the absence of mixing, 
x' = = and R(tlT) = Rd- 

Our measurement uses data collected by the CDP II 
detector at the Permilab Tevatron collider, from February 
2002 to January 2007, corresponding to an integrated lu- 
minosity oi~ 1.5 fb~^ forpp collisions at ^/s = 1.96 TeV. 
CDP II [l^ is a multi-purpose detector with a magnetic 
spectrometer surrounded by a calorimeter and a muon 
detector. The detector components pertinent to this 
analysis are the silicon microstrip vertex detector, the 
multi-wire drift chamber (COT), and the 1.4 T magnet 
which together measure the trajectories and momenta of 
charged particles. The COT measures ionization energy 
loss for a charged particle, which is used for particle iden- 
tification (PID). 

Events are selected in real time with a trigger system 
developed for a broad class of heavy ffavor decays. The 
trigger requirements used here are the same as those 
described for our previous measurement of the time- 
integrated value of R [3], which used a smaller data 
sample. The trigger selects events with a pair of oppo- 
sitely charged particles that are consistent with originat- 
ing from a secondary decay vertex separated from the 
beamline. 

In the off-line analysis, we reconstruct the "right-sign" 
(RS) CP decay chain D*+ tt+D", D° K-Tr+, 
and the "wrong-sign" (WS) decay chain D*+ tt+D", 
D'^ K^TT^. The relative charges of the pions deter- 
mine whether the decay chain is RS (like charge) or WS 
(opposite charge). The reconstruction method is similar 
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to that used for our previous time-independent measure- 
ment. The RS and WS D* decays have the same kine- 
matics, and differ only in decay time distributions. To 
reduce systematic uncertainties, we use the same selec- 
tion criteria (cuts) for both the RS and WS decay modes. 
Analysis cuts were optimized before the WS candidates 
were revealed, and were chosen to maximize the expected 
WS signal significance. 

The candidate reconstruction starts with a pair 
of tracks from oppositely charged particles that satisfy 
the trigger requirements. The tracks are considered with 
both K~TT~^ and tt^K^ interpretations. A third "tag- 
ging" track, required to have pr > 0.3 GeV/c, is used 
to form a D* candidate when considered as a pion and 
combined with the candidate. 

We apply two cuts to reduce the background to the 
WS signal from RS decays where the decay tracks 
are misidentified because the kaon and pion assignments 
are mistakenly interchanged. As determined from the 
data, 96.4% of decays with correct mass assignment 
are reconstructed with Kn invariant mass mxn- within 
20 MeV/c^ of the mass. The rriKTr distribution for 
misidentified decays is much broader, and has only 
22% of the events within the same mass range. We re- 
move WS candidates that have a RS mass within that 
range. This cut excludes 96.4% of RS decays and retains 
78% of the WS signal. We also impose a cut based on 
PID which is used to distinguish pions from kaons for all 
three tracks in the decay chain. This cut, described in 
Ref. [3], further helps to reject misidentified decays. 

We use a series of cuts based on the decay topology 
of signal events in which a D* is produced at the pri- 
mary vertex, the tagging pion also originates from the 
primary vertex, and the travels a measurable distance 
before decay. The vertex-based cuts reduce combinatoric 
background from improper combinations involving one 
or more tracks that do not originate from the D* de- 
cay chain of interest. We require the transverse decay 
length significance Lxy/dxy to be greater than 4, where 
L'xy — ^-pt/pt, r \s the distance between the primary 
and decay vertices, pt is the transverse component 
of the momentum of the candidate with respect to 
the beamline, and Uxy is the uncertainty on L^y. The 
tagging pion track must have do < 500 /xm, where the 
transverse impact parameter do is the distance of clos- 
est approach between a track and the primary vertex in 
the plane transverse to the beamline. The tagging pion 
must also have a point of closest approach to the primary 
vertex less than 1.5 cm along the beamline. 

The ratio t/r is determined for each D'^ candidate by 
t/r ~ ra^ioLxy/ [ptt), where mjja = 1.8648 GeV/c^ and 
T — 410.1 fs are the world average values for the in- 
variant mass and lifetime, respectively [lO]. To study 
R{t/T), we divide the data into 20 bins of t/r rang- 
ing from 0.75 to 10.0, choosing bins of increasing size 
from 0.25 to 2.0 to reduce statistical uncertainty at larger 
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FIG. 1: Time-integrated distribution for "wrong-sign" D 
K^-K~ signal yield as a function of Am. Also shown is the 
result of a least-squares fit using an empirical function for the 
signal (dark shaded region) and a power law for the back- 
ground (light shaded region). 



times. The bin sizes are larger than the t/r resolution of 
w 0.16. 

After RS and WS candidates are separately divided 
into t/r bins, they are further divided into bins of mass 
difference Am = mK-mr — ttiktt — . For each Am bin, 
we perform a binned maximum likelihood fit of the cor- 
responding niK-n distribution to determine the signal 
yield. The distribution of D'^ signal yield versus Am is 
fit using a least-squares method to get the D* signal for 
each time bin. The D* fit procedure is illustrated by the 
time-integrated WS Am distribution shown in Fig. [TJ 

The signal shapes for the individual m^^ and Am 
distributions are fixed from the RS time-integrated fits. 
For each rriK-n distribution, a parabola with fioating pa- 
rameters is used to fit the background. The background 
shapes for all the Am WS (RS) distributions are fixed to 
the shape determined for the time-integrated WS (RS) 
distribution. The amplitudes of the signal and back- 
ground shapes are determined independently for all rriK-n 
and Am fits. The RS distributions have similar amounts 
of background as the WS distributions, but the RS signal 
is about 250 times larger. 

The D* mesons that originate from beauty hadron {B) 
decays must be treated as background to avoid the com- 
plication of measuring the decay length from the 
B decay point instead of the primary vertex. The D* 
mesons produced promptly at the primary vertex have 
a narrow do distribution, with a shape independent of 
t/r. The background from non-prompt D* mesons from 
the decay chain B D* D'^ have a broad do distri- 
bution, due to the decay length of the B hadrons. The 
width of the broad distribution increases with increasing 
t/r. An example do distribution is shown in Fig. [2l The 
shapes of the prompt and broad distributions are deter- 
mined from RS data. The WS shapes are the same as 
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FIG. 2: The distribution of transverse impact parameter do 
for 5 < t/r < 6 for "right sign" background-subtracted D* 
mesons. The result of a binned maximum likehhood fit shows 
the narrow peak due to promptly produced D* mesons (dark 
shaded) and the broad distribution due to non-prompt D* 
mesons from B decay (hght shaded). 
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FIG. 3: Ratio of prompt D* "wrong-sign" to "right sign" 
decays as a function of normalized proper decay time. The 
dashed curve is from a least-squares parabolic fit, which de- 
termines the parameters Rn, y' , and a;'^. The dotted line is 
the fit assuming no mixing. 



the RS shapes. For each of the 20 t/r bins, the prompt 
WS (RS) signal is determined from the number of WS 
(RS) D* mesons and the shapes of the do distributions. 
The ratio of non-prompt to prompt signal is ~ 0.02 at 
t/r ~ 2 and increases with increasing t/r due to the 
faster exponential fall-off with t/r for compared to 
B. At t/r = 7, the ratio is « 1. 

The time-integrated prompt D* signals are (12.7 ± 
0.3) X 10^ WS events and (3.044± 0.002) x 10^ RS events. 
The ratios of prompt WS to RS signal for the 20 t/r 
bins are shown in Fig. [3] The uncertainties for each 
bin include statistical and systematic contributions. The 
significant systematic uncertainties are due to the back- 
ground shapes for the mK-n, Am, and c?q distributions, 
which are described by parameters that are allowed to 
vary in the fitting procedure. We used simulation to 
confirm that our choice of decay time bins does not sys- 
tematically affect the result. The detector acceptances 
for RS and WS decays are nearly identical, and their dif- 
ference contributes a negligible systematic uncertainty in 
the ratio R. The large uncertainty in the first t/r bin is 
due to low signal statistics caused by the trigger turn-on 
at small t/r. In the last two bins, the uncertainties are 
large because the exponential fall-off of the WS signal 
with t/r results in smaller numbers of signal events. 

A least-squares parabolic fit of the data in Fig. [3] to 
Eq. Il]) determines the values and uncertainties for the 
parameters Rd, y', and a;'^, which are listed in Table [D 
Since the value of x'^ is unphysical (less than zero) , but 
consistent with zero, we also fit the data with the con- 
straint x'^ — 0. The values of Rd and y' are consistent 
with and without the constraint. The values and preci- 
sion of the parameters measured by CDF are comparable 
to those from the best previous measurements, as shown 



TABLE I: Fit results for the R{t/r) distribution. The uncer- 
tainties include statistical and systematic components. The 
correlation coefficient between y' and x'^ for the uncon- 
strained fit is —0.98. The no-mixin g fi t is consistent with 
our previous time- independent result [19l |. 



Fit type RpjlO-'') y' jW') a:'^ (10~^) d.o.fT 

Unconstrained 3.04 ± 0.55 8.5 ± 7.6 -0.12 ± 0.35 19.2 / 17 
Physically 

allowed 3.22 ± 0.23 6.0 ±1.4 19.3 / 18 

No mixing 4.15 ±0.10 36.8 / 19 



in Table El 

To determine the consistency of our data with the no- 
mixing hypothesis, we compute Bayesian contours con- 
taining the region with the highest posterior probability. 
The probability density is calculated as the product of 
a likelihood C and a prior, divided by a normalization 
factor. The likelihood is £ = exp(— x^/2), where is 



TABLE II: Comparison of the CDF result with recent mea- 
surements. All results use —^ K^-k~ decays and fits as- 
suming no CP violation. The uncertainties include statistical 
and systematic components. The significance for no mixing is 
given in terms of the equivalent number of Gaussian standard 
deviations. 



Mixing 

Experiment fic(10~") y' [IQ'^) a:'^ (10~^) Signif. 
CDF 3.04 ± 0.55 8.5 ± 7.6 -0.12 ± 0.35 3.8 

BABAR [8] 3.03 ± 0.19 9.7 ± 5.4 -0.22 ± 0.37 3.9 
Belle [9] 3.64 ± 0.17 0.6 tfl 0.18 +023 2.0 
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FIG. 4: Bayesian probability contours in the x''^ — y' pa- 
rameter space corresponding to one through four equivalent 
Gaussian standard deviations. The closed circle shows the 
unconstrained fit values for the mixing parameters. The open 
diamond shows the values from the physically allowed fit 
{x'^ > 0). The cross shows the no-mixing point. 



are needed to explain this phenomenon. 
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computed from the data in Fig. [3] for a particular set of 
fit parameters. The normalization factor is the integral 
of C over the parameter space. A flat prior is used for all 
three parameters, and is treated as a Bayesian nui- 
sance parameter. The contours are insensitive to modest 
changes in the prior. The contours in the x''^-y' plane 
are shown in Fig. (4] The no-mixing point lies on the 
contour which excludes a region containing a probability 
of 1.5 X 10~^, equivalent to 3.8 Gaussian standard devi- 
ations. We also computed contours with the constraint 
x'^ > and find a probability for no-mixing consistent 
with the value obtained without the constraint. 

We tried alternate procedures to determine the prob- 
ability for no mixing. We fit the data in Fig. [3] with the 
constraint y' = x'^ = 0, with results as given in Table [J 
The change in log likelihood (— 2Aln£) between the un- 
constrained and no-mixing fits has an approximately chi- 
square distribution for two degrees of freedom. From 
Table m — 2Aln£ = 17.6, which corresponds to a prob- 
ability of 1.6 X lO^'*. We also made a frequentist check 
using ensembles of simulated Rit/r) measurements with- 
out mixing. The probability for a simulation to have a 
value of -2Aln£ > 17.6 is 1.3 x lO"'^. The probabil- 
ities from both of these checks are consistent with that 
obtained using Bayesian contours. 

In conclusion, our data shows evidence for — l)" 
mixing in the K^n" channel, providing the first con- 
firmation of the evidence from the BABAR experiment. 
The mixing could be due to standard model long-range 
intermediate states or due to new physics. Improved re- 
liability of standard model calculations and future mea- 
surements of mixing signatures with improved precision 
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